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We disprove the Vafa-Witten theorem on the impossibility of spontaneously breaking parity in 
vector-like gauge field theories, identifying a mechanism driven by quantum fluctuations. With the 
introduction of a meromorphic Lattice formulation, deflned over 5 dimensions, we demonstrate that 
the minima of the free energy can be distinct from the maxima of the partition function : identifying 
and evaluating a suitable contour for the partition function defined such that asymptotic behaviour 
of the complex action is non-oscillatory. 
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Recently attempts have been made to determine the 
properties of the Grand Canonical Partition Function of 
QCD at nonzero baryon densities using numerical Lat- 
tice QCD approaches 0. These results are directly rele- 
vant for determining the equation of state of QCD for ex- 
perimental Heavy Ion Collider results. These efforts arc 
frustrated numerically by the difficulties in constructing 
efficient sampling algorithms at finite baryon densities. 
One faces the problem of a complex-valued Lattice ac- 
tion, whose phase oscillates rapidly as a function of the 
configuration within a Lattice statistical ensemble. This 
then prevents the effective use of importance sampling 
methods. All recent sampling improvement treatments in 
this vein appear now to suffer from the same generic flaw 
of an unphysical light pseudo-baryonic transition |3, Q • 

To gain understanding of this feature an analytic treat- 
ment of partition function zeroes for QCD and QCD-like 
theories with a complex action problem was recently pro- 
posed 0, 13 . Here the properties of the zeroes are iden- 
tified as relating directly to specific regions of parame- 
ter space in the continuum where the partition function 
is itself ill-defined. This is linked in the discussion to 
the Vafa-Witten Theorem on the impossibility of spon- 
taneously breaking parity symmetry in QCD 0. A lo- 
cal bosonic parity breaking source term in the action is 
shown to generate a similarly pathologically ill-defined 
continuum partition function. This is demonstrated ex- 
plicitly through the properties of the zeroes, and is used 
to justify the main theorem result. Although not in- 
cluded in this specific zeroes analysis the Vafa-Witten 
theorem result is understood to represent a wider pathol- 
ogy in Lattice simulation 0, Hi . Bosonic terms with sub- 
structure in the Vafa-Witten formalism give rise to un- 
physical diagrammatic contributions to the Lattice QCD 
partition function. For example, in Nf = 2 QCD with a 
four fermion operator source , or the Aoki phase with 
Wilson fermions [Hill El- 
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We now apply a new quantum system partition func- 
tion zeroes treatment to this problem, and show that in 
contrast it is possible to give a well-defined continuum 
partition function for all generic complex action prob- 
lems (of those of the form in the Vafa-Witten discus- 
sion) for all parameter space regions. Under the proviso, 
that, the local source term has a highly nonlocal com- 
ponent for which the continuum limit has been reached 
over an additional phase space dimension. We introduce 
a new quantum loop operator definition for the partition 
function and show that it is formally equivalent to the 
Vafa-Witten construction, but in addition, meromorphic. 
This new analysis is used to rigorously disprove the Vafa- 
Witten result, and demonstrate that it is in fact possible 
for any global Lattice symmetry (like taste) to be sponta- 
neously broken via the system dynamics. Specifically, we 
demonstrate that the branch cut argument introduced 
in 0, 1^ ^3 ^3 used to argue that the partition func- 
tion is ill-defined is specious. The branch in fact can 
be resolved analytically and unambiguously for all Lat- 
tice system sizes, N , and in particular in the asymptotic 
limit. 

The new analysis also indicates that unless the Lat- 
tice simulation of a quantum system with a complex ac- 
tion problem produces extremely finely balanced statis- 
tical averaging between degenerate vacua (namely that 
which in the asymptotic limit N oo connects smoothly 
with the continuum) then spurious phase behaviour will 
be observed. We have applied a new fine-balancing 
scheme recently, based on this understanding, via numer- 
ical simulation to several Lattice QCD models at finite 
baryon densities and independently quantum spin sys- 
te ms ( which have a hidden continuous gauge symmetry) 
H HH EIHIIl- This article then, in part, provides 
a generic proof of the existence of the continuum limit 
of these ongoing Lattice analyses, and the validity of our 
generic continuum matching program for these calcula- 
tions. 

The purpose of this new article is to establish the origin 
of large light disconnected contribution in recent lattice 
QCD schemes at finite baryon densities in a generic man- 
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ner. Specifically, we wish to establish the relationship 
between the statistical ambiguities of numerical sampling 
and the parity symmetry violation arguments of Vafa and 
Witten. This is expressed through a saddle point point 
discussion of the partition function and is specific to our 
new quantum loop operator formulation. To do this, fol- 
lowing Lee and Yang , we wish to now argue for sep- 
arate treatment for the approach to the continuum limit 
of the zeroes polynomial in quantum systems, as opposed 
to classical systems. We therefore rederive the partition 
function zeroes polynomial proof of Lee and Yang and 
the validity of its approach to the continuum limit for the 
quantum case, formulated in terms of loop operators. In 
doing so, we demonstrate that this is the generic origin 
of these disconnected contributions in numerical simu- 
lations. This then has direct applicability in the design 
of new improved algorithms - to include correct quantum 
effects - both for Lattice QCD and quantum spin systems 
that are otherwise spoilt by multiplicative nonrenormal- 
izable anomalies induced through the noncommutativity 
of finite Lattice system dynamics [23, ■ 

This article is then organised as follows : in Section 
^we introduce a new quantum loop operator formalism 
for the partition function and identify that this is for- 
mally equivalent to the local bosonic operator expansion 
on finite system size introduced in the Vafa- Witten dis- 
cussion. We then consider the continuum limit of this 
formulation and demonstrate that the zeroes of the loop 
operator polynomial are formally equivalent to partition 
function zeroes, but only as defined with reference to a 
prior noncommuting continuum limit being taken over 
an additional dimension in an extended phase space. We 
then derive the saddle point solution for the partition 
function of this complex-valued quantum loop operator 
formulation. 

In Section m we discuss the analytic continuation of 
the saddle point solution with reference to the symme- 
try breaking properties of the newly defined quantum 
partition function, and the Vafa- Witten discussion. We 
demonstrate that the results can be analytically contin- 
ued across the entire phase plane in the zeroes expansion 
parameter. Therefore any nonperturbative operator sat- 
isfying the quantum loop operator definition can be ana- 
lytically continued throughout the complex phase-plane, 
and is an exact analytic solution to the complex action 
problem. 



I. QUANTUM LOOP OPERATORS 

The partition function, for a general continuous 
group G we define in terms of quantum loop operators 
by, 



Z]^(t)= J dg exp J Ag{ns) z - V{ns) ds 
with t e M, z e C, 



(1) 



N(g)T 



(n|n) 



where Xija{n), X'^j^ln) E C, {a} form some finite sub- 
set of the elements of G, and g is a general element of 
G. We treat the complex-valued continuous symmetry 
generalisation of the problem considered in j^i 2^ ■ The 
fact that the operator definition is diagonalisable implic- 
itly defines a local gauge transformation on the Hilbert 
space, 



N<g)T 



E E ^^3o{ri) zTTtH Is) 



= (n'|.9)(n'|^|n)(n|5)-^ = (n'|A,|n>. (4) 

We should compare this directly with the single plaque- 
ttc model given in terms of the Vandermondc determi- 
nant discussed in |23 . l25l |. which yields a similar form 
of time-ordered product. The connection to Vafa- Witten 
formalism is made explicit by noticing that equivalently 
the integration of the exponent in can be performed over 
the spatial Lattice volume via a change of basis, n'^. The 
partition function is then given in the usual Vafa- Witten 
form, 



y dg exp y Agiris')^ ~y[ns')ds' 



Zt{z) = I dg exp 
= / dg exp 



A' {n,,)it~V'{n,,)ds' 



(5) 



A'gM = lm[AgM], (6) 
V'M = VM-Rc[Agins')], (7) 

where V{ns') E M. since it is hcrmitcan in the previous 
basis. Us- 



A. Thermodynamic Limit Monotonicity 

It is surprisingly simple to prove the existence of the 
thermodynamic limit for partition function in The 
proof consists of two parts. The first step is to identify 
the saddle point solution of the partition function. The 
second step is then to determine the limiting value of 
this asymptotic solution when the neighbourhood of the 
saddle point is extended to cover the entire phase space. 
Apart from isolated singularities the partition function is 
then necessarily holomor phic. T his is the essence of Lee 
and Yang's discussion in 
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For the quantum case we face now an additional sub- 
tlety. This has not been previously addressed in the lit- 
erature. Partition functions for quantum systems are de- 
fined over an extended phase space - a trace is taken 
over an additional dimension. Two examples are given 
above in and 101 ■ For a quantum system to be truly 
holomorphic, and the thermodynamic limit exist, it is 
therefore not merely sufficient that the thermodynamic 
limit be proven to exist over either ris or n^'. Rather, 
the thermodynamic limit must exist over n. Expressing 
this another way, existing partition function studies in- 
voking this holomorphic property over ris' (holomorphic 
in z) have a thermodynamic limit that is strictly speak- 
ing only defined up to infrared divergences (ie. they are 
not necessarily simultaneously holomorhpic in t) . For Lee 
and Yang with classical systems, which have no trace re- 
lation, this complication docs not arise. The partition 
function in this case is given directly over n, without the 
time-ordering of the integrand. 

We now prove the existence of the thermodynamic 
limit for the new quantum loop operator formulation for 
the partition function in JQ). The saddle point, which 
maximises the partition function in some region of the 
extended phase space M, is defined through. 



Rm = sup i?(m) . 

M 



(8) 



where m is parameterises some contour m N and R 
is the integrand of the partition function. The value of 
the integrand in the neighbourhood of the saddle point, 
Ms, is defined through. 



Rsim) = R{m) - Rm < . 



(9) 



Not every point contained in A/ is necessarily in the 
neighbourhood of the saddle point. We consider how 
the neighbourhood evolves as the contour is increased in 
length : the limit where m ^ n. 

lim Rs{m)^ lim R{m) - Rm < ■ (10) 

m — >n m — >n 

The relation between the consecutive basis elements 
(along m) is given partly through the local gauge sym- 
metry transformation in However, we can have a 
situation for some general group G in which we have 
multiple degenerate solutions of ifTUI) . In this case a sub- 
space must be projected out to time-order the integrand 
and resolve this degeneracy. With this resolution is then 
necessarily monotonic and decreasing, since the potential 
nonlocal ambiguities are removed. Next, taking the limit 
liniA/^AfigjT we have, 



lim [ lim Rs(m)] 

lim [ lim R(m) - Rm] < ■ (H) 

This limit is also monotonic and decreasing with above 
choices, and therefore the thermodynamic limit of the 



quantum loop operator partition function in is well- 
defined. The holomorphic properties we are therefore 
interested in are those of the space ® T, modulo, the 
time-ordering necessary to resolve the nonlocal ambigu- 
ities of the gauge group. It should be clear now from a 
numerical perspective that it is not necessary to perform 
either integral in the integrands exactly : in order to be 
simultaneously holomorphic in t and z it is sufficient for 
the integrand to cover N ®T. 



B. Saddle Point Determination 

We are interested to find an explicit expression for 
the evolution of the partition function, Zn, in the limit 
N oo, N > 0. There are exactly N singularities in 
the partition function in as a consequence of the above 
monotonicity arguments. By construction, the loop op- 
erator definition of Zf^ is meromorphic so we can equally 
well formulate the discussion in terms of the poles of Zm, 
{Vfe}, or the zeroes, {Xk}- The poles turn out to be eas- 
ier to treat than the zeroes, but for clarity the relation 
between the poles and zeroes is given through. 



1 

2TTi 



A{z 
A'{z 



N 

n 

fc=0 



z + Xk 
z + Vk 



(12) 



where n and n' are the respective homotopies of the ze- 
roes and poles, and C is some contour homologous to 
in TV (g) r. This property will become important shortly 
when we consider finite multiplicities for the singularities, 
or equivalently, global symmetries in N ^T. 

The schema we now follow is the complex- valued gen- 
eralisation of the problem considered in [23, 123 . Namely, 
we define the Laplace transform of the partition function, 
Zn, as a function of the poles, {Vfe}. Then we formulate 
the inverse Laplace transform that gives us back the par- 
tition function expressed as a function of t, N and {Vfc}. 
Finally, we evaluate this integral equation for Zpf via 
standard asymptotic series approaches in order to iden- 
tify the dominant behaviour of the partition function in 
the above limit. In particular, this approach allows us to 
identify the saddle point solution of the partition function 
as distinct from the asymptotic convergence properties of 
the zeroes towards the Re z axis. 

The Laplace transform of the partition function poly- 
nomial in at finite is given in terms of the poles of the 
partition function by. 



^ poo ^1 

Z^iz) = / dte-^'Zr,{t) = e^l[ — rr, (14) 
Jo k=o + ^'^ 

where z, e, Vk G C. The inverse Laplace transform of is 
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given by, 



ZN{t) 



dz exp [Nip(z)] , 

'c 

(15) 



generic form above in. Solving explicitly for (p'{zo) 
we have, 



N 

E 

k=0 



zo 



Vk 



t. 



(17) 



ip{z) 



zt 
N 



^ 1 

A:=0 



z + Vk 



(16) 



The contour C in, is defined as any suitable path running 
between ±ioo, for which any singular behaviour is appro- 
priately treated. We want to now evaluate the asymp- 
totic behaviour of the partition function separately for 
the two cases : when a parity symmetry is, and is not, 
realised on the vacuum. 

In the first case we have no special symmetry relations 
defined between the poles, and (p{z) is therefore of the 



Note, the location of the saddle point, zq, is a function of 
both {Vk} and t. Potentially the saddle point lies to the 
left of some subset of the poles, {T^m}, m < N. Therefore 
the integral in is given by, 

= — / exp[N^{z)]dz + y] — Yi 



1 



y- 



(18) 

where Cz^ is the contour that passes through zq- For 
this first case, when one assumes no parity symmetry is 
realised with the gauge fields, the asymptotic behaviour 
of the partition function is therefore given by. 



hm ZN{t) = — — 



exp 



exp 



Nip{z^) ~ -9{zo) 



-2tt 



N\ip"izo)\ 



N 

zot - ^ log 

k=0 



Zo + Vk 



fizo) 



+ C»(7V-3/2) 



2TTe^N\ip"{zo)\ 



0{N 



-3/2N 



(19) 
(20) 



where 9{zo) = Argip"{zo). 

For the second case it is straightforward to define a 
parity symmetry action on Zt^{t) simply by reversing the 
orientation of the elementary plaquettes in N ®T . This 
then yields the symmetry of the zeroes {Afc,A^.}. From 
the Argument Principle in (|13|l the poles are then given 
through. 



Zn{z) 



N 

n 



\^{z + Vk){z + V*) 



Lp{z) 



f 1 / 



^{z + Vk){z + V^) 



(21) 



(22) 



The singularities at z = Re ± Im 14 therefore cor- 
respond to branch points, with a branching exponent of 
— 1/2. Again, we want to define the contour such that it 
passes through zq, but this is again potentially obstructed 
by singularities lying to the right of zq. We therefore 
choose the branch cut to correspond to this deformation 
over the saddle point. If Im Vfc < the deformation will 
start at the branch point and descend to —ioo. If the 
Im Vk > the deformation will start at the branch point 
and ascend to zoo. If Im = the singularity is a pole, 
which we will treat as a residue contribution as before. 



This choice of branch cut does not introduce an arbi- 
trariness into the calculation. Although it is not now 
possible to analytically continue across the cut, ip{z) is 
specified unambiguously along the original contour, and 
the analytic continuation between contours is also speci- 
fied unambiguously. Crucially we keep the same branch 
cut for the pairs {VJ, V;*}. The cuts to the left and right 
of the contour deformation C are defined through. 



(z + Vi)-'/^ = e--'\z + Vi)Z'/\ 



(23) 



The analogue of the residue contribution for the branch 
points is therefore. 



gJV^(.,MO(^ + Vn-'^^iz + Vi)-_^'^dz 
gJV^(z,MO(^ + VCr^'\z + Vi)-^'^dz 



(24) 

2 j ^N^(z,k^i)^^ _^ VC)-^'^{z + Vi)Z^'^dz 

(25) 
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where, 



b{0)^{-2ilmVi.^'{Vi,k^l)) (27) 



1 



Re(t^z + Vfe) 



(28) 

and the 6"(0) form the analogue of the Laurent series for 
the branch points. The solutions are not oscillatory in 
TV, as was found with previous discussions. In our new 
formulation we can see as well that the behaviour of the 
branch points can wholly dominate the asymptotic form 
of the partition function in the case in which a parity 
symmetry is realised, solely dependent on the relative 
magnitude of the prefactor, cxp[tp(yk-,t)]. 

We are able to reach this new point of unambiguous 
saddle point solutions for two reasons. Firstly, the new 
problem is defined over a higher dimension than previ- 
ous discussions, so z and t are separated. This means 
that the singularities can be distinguished from asymp- 
totic scaling. Secondly, because the new construction is 
meromorhpic we can identify all the singularities in Zn, 
which are of a finite number. 

The above arguments also are readily generalised to 
cases with different branching exponents - corresponding 
to different global symmetries. Specifically, if the branch- 
ing exponent is less than -1 then the above asymptotic 
identifications for the branches can be performed via in- 
tegration by parts. 



II. SPONTANEOUS SYMMETRY BREAKING 

We have constructed in a partition function that 
is formally equivalent to that of the Vafa-Witten discus- 
sion. It consists of a local bosonic path integral, on which 
is defined a symmetry breaking source term, Aglngi). 
To relate our saddle point identification further to the 
Vafa-Witten discussion on spontaneous parity violation 
we now need to evaluate the vacuum expectation of the 
source term. Since the path integral is time-ordered this 
is given by. 



(O) = lim 



1 d log {Zn it)) 



N->oo \ N 



dt 



lim 



I Jdz ze^* Zn{z) 



(29) 



(30) 



N^oo \yN Jdz e^*ZAr(z) / 

We therefore have two possible solutions, from H22|l and 
(|26|l . for the case when parity is realised by the system. 



In particular we can see for this second case in (|26|l . asso- 
ciated with the branch point, that the free energy density 
minima does not correspond with the partition function 
maxima at zq. What we can see through our new for- 
malism is that at finite N there is a finite subtraction 
to made to properly define the contour for the partition 
function, and this itself can give rise to a parity breaking 
contribution. From H22|l at most one of these expecta- 
tions can be zero simultaneously, and so the potential for 
spontaneous symmetry breaking is well-defined. 

Our new analysis does not mean that it is particularly 
likely that parity be spontaneously broken in QCD, since 
in our formulation it is identified as a quantum effect 
driven by the system dynamics through {14}. Plausibly, 
these symmetry breaking effects are only observed in in- 
stances where the statistical properties of the the Lattice 
ensemble become very much non quasi-classical, ie. the 
orders of the contributions in (|22|l and H26|l become com- 
parable lUllS- 

As a second result we should notice from that it is 
possible to analytically continue the partition function 
throughout the entire z plane. By definition the partition 
function is holomorphic in this region. We can deduce in 
principle it is therefore possible to resolve the complex 
action problem via numerical simulation. If a lattice de- 
termination correctly determines the polynomial expan- 
sion in up to quantum fluctuations at any single point in 
z, ie. that generates a covering of the Hilbert space in 
^ and 121), the system can be smoothly mapped to all 
point in the z plane. Since the result smoothly maps to 
the continuum the Lattice evaluation is also sufficiently 
defined for the dynamics to generate a nonzero vacuum 
expectation value of the source, O, and to observe spon- 
taneous symmetry breaking in the continuum. 

Conversely, we can understand the origin of the discon- 
nected contributions in existing Lattice evaluations. Al- 
gorithms can be based around the principle that a global 
symmetry is exactly held, and equal Monte Carlo weights 
assigned for each branch. However, because only the pos- 
itive definite contributions are kept, corresponding to real 
probabilistic weights, the opposite cut is in effect taken 
for each branch point multiplet in (|21|l . Genuine quan- 
tum effects can never therefore be observed as the can- 
cellation between sectors is exact configuration by con- 
figuration, by construction and furthermore the system 
is nonrenormalizable l3lLl33l. 



{O)=zo or {0)=Vi. 



(31) 



III. SUMMARY 

We have introduced a quantum loop operator formal- 
ism for the partition function of vector-like gauge field 
theories, and identified that this is formally equivalent 
to the local bosonic operator expansion on finite system 
size introduced in the Vafa-Witten discussion. The Lat- 
tice formulation is given over an extended phase space 
N ® T, and the source parameter is itself discretized in 
the new approach. We have then considered the contin- 
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uum limit of this space and identified that the zeroes of 
the loop operator polynomial are formally equivalent to 
partition function zeroes. Although this result is distinct 
from classical systems requiring a prior noncommuting 
limit to be taken over T to ensure the formulation is 
meromorphic. The saddle point solution for the asymp- 
totic behaviour of this complex-valued partition function 
was then determined. Spontaneous symmetry breaking 
is well-defined for case in which a parity symmetry is re- 
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